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We investigated the origin of perpendicular magneto-crystalline anisotropy (MCA) in Llo ordered 
FeNi alloy using first-principles density-functional calculations. We found that the perpendicular 
MCA of Llo-FeNi arises predominantly from the constituent Fe atoms, which is consistent with 
recent measurements of the anisotropy of the Fe orbital magnetic moment of Llo-FeNi by x-ray 
magnetic circular dichroism. Analysis of the second-order perturbation of the spin-orbit interaction 
indicates that spin-flip excitations between the occupied majority-spin and unoccupied minority-spin 
bands make a considerable contribution to the perpendicular MCA as does the spin-conservation 
term in the minority-spin bands. Furthermore, the MCA energy increases as the in-plane lattice 
parameter decreases (increasing the axial ratio c/a). The increase in the MCA energy can be 
attributed to further enhancement of the spin-flip term due to modulation of the Fe d(xy) and 
d(x 2 — y 2 ) orbital components around the Fermi level under the compressive in- plane distortion. 
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I. INTRODUCTION 

Ferromagnetic materials with strong magneto- 
crystalline anisotropy (MCA) have attracted much 
attention for application in high density magnetic 
recording media and nonvolatile magnctorcsistivc 
random access memory (MRAM). In order to ensure 
sufficient endurance against thermal fluctuations when 
downsizing a recording bit or memory cell, it is inevitable 
to adopt ferromagnetic materials with higher MCA, 
such as Llo-ordered FePtji CoPt^ and ZX)22-ordered 
Mn3_aGai^ Furthermore, perpendicularly magne- 
tized electrodes of magnetic tunnel junctions (MTJs) 
lower the critical current density for current-induced 
magnetization switching by spin-transfer torque^ 

Recently, Llo-ordered FeNi, which is free from rare- 
earth and/or noble metal elements, have been suc- 
cessfully fabricated as thin films using an alternative 
monatomic layer deposition technique.— ~— The MCA en- 
ergy of Llo-FeNi thin films depends on the sort of buffer 
layers 7 - as well as the degree of chemical order 4 Fur- 
thermore, x-ray magnetic circular dichroism (XMCD) 
and magneto-optical Kerr effect measurements were car- 
ried out on alternately layered FeNi thin films grown 
on Ni/Cu(001) substrates^ By analyzing the Fe L-edge 
XMCD spectra using the sum rule ; 10 ' 11 they obtained the 
anisotropy of the orbital magnetic moment, and found 
that the Ni-sandwichcd Fe layer has a MCA energy of 
10 /xeV, while the Fe-sandwiched Ni layer has a positive 
MCA energy of 60 ^cV. 

The MCA of ordered Ll -type alloys TX (T = Mn, 
Fe, Co, Ni, X = Au, Pt, Pd, Co, Ni) have been theoreti- 
cally well studied using first-principles density functional 
calculations.— ~— For Llo-type FePt or CoPt, the strong 
perpendicular MCA was attributed to the Pt atom due 
to the strong spin-orbit interaction, in which Fe or Co 



atoms only induce the spin magnetic moment on the Pt 
siteJ^ On the other hand, Llo-FeNi also shows a per- 
pendicular MCA, despite the weak spin-orbit interaction 
of Ni compared to that of Pt. Therefore, the physical 
origin of the perpendicular MCA of Llo-FeNi remains 
unclear at present. Furthermore, recent XMCD mea- 
surements reported the anisotropy of the orbital mag- 
netic moment These reports discussed the MCA on 
the basis of Bruno's relation^ in which the MCA en- 
ergy is proportional to the difference in orbital magnetic 
moment between two magnetization directions. How- 
ever, there are additional terms related to spin- flip exci- 
tations between the exchange-splitting majority-spin and 
minority-spin bands in the second-order perturbation of 
the spin-orbit interaction. The effects of the spin-flip 
terms on the MCA energies can be neglected for systems 
with strong exchange coupling, such as a free-standing Fe 
or Co monolayer^ while the spin-flip term makes a large 
contribution to MCA for systems including nonmagnetic 
elements. In fact, previous theoretical work revealed an 
enhancement of the Pt orbital magnetic moment of Llrj- 
FePt in the in-plane magnetization direction; 12 ' 1 while 
real-space analysis of the MCA energy indicated a strong 
perpendicular contribution from Pt atoms This means 
that the spin-flip term is essential to understanding the 
perpendicular MCA of ordered Llo-type alloys. 

In this paper, we discuss the MCA energy of Llo- 
FcNi theoretically on the basis of first-principles density- 
functional calculations in order to elucidate the physical 
origin of perpendicular MCA in Llo-FeNi. In particu- 
lar, we focus our attention on the spin and atomic-site 
dependence of the MCA energy using a second-order per- 
turbation of the total energy due to the spin-orbit inter- 
action. Furthermore, we discuss the dependence of the 
MCA energy on the in-plane lattice constant of Llo-FeNi, 
providing useful guidelines for the choice of buffer layers 
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In order to obtain a fundamental understanding of the 
origin of MCA in Llo-FeNi, we consider the second-order 
perturbation of the total energy due to the spin-orbit 
interaction i?sO ) 21l32 ' S3 which is given by, 
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FIG. 1: (color online) Schematic figures of the atomic arrange- 
ment in Llo-FeNi. The coordinate system and the primitive 
tetragonal unit cell used in the present study are represented 
by bold solid lines. 



in Llo-FeNi thin films to realize larger uniaxial MCA 
constants K u above 1.0 MJ/m 3 . 



where \kna) is an unperturbed state of energy ej£L with 
indices of fe-point k, band n, and spin a. In the tight- 
binding regime for spin-orbit coupling, Hgo is given by 
the sum of the contributions from each atomic site 
due to the localized character of the spin-orbit coupling 
constant £j, 
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II. COMPUTATIONAL PROCEDURE 

A. First-principles calculations 

We carried out first-principles density-functional cal- 
culations using the Vienna ab initio simulation pack- 
age (VASP)i 24 ' 25 For the exchange and correlation en- 
ergy, we adopted the spin-polarized generalized gradi- 
ent approximation (GGA) proposed by Perdew, Beckc, 
and Ernzcrhof (PBE)^ The nuclei and core electrons 
were described by the projector augmented wave (PAW) 
potentia l 27 ' 28 and the wave functions of valence electrons 
were expanded in a plane-wave basis set with a cutoff 
energy of 336.9 eV. The fc-point integration was per- 
formed using a modified tetrahedron method with Blochl 
corrections 2 ^ with 19800 fc-points in the first Brillouin 
zone corresponding to the primitive unit cell of Llo-FeNi. 
The spin-orbit interaction was included through the force 
theorem^ We determined the MCA energy from the dif- 
ference in the sum of energy eigenvalues for magnetiza- 
tions oriented along the in-plane [100] and out-of-plane 
[001] directions. The MCA energy was defined to be pos- 
itive when the magnetic easy axis was perpendicular to 
the plane. As the unit cell of Llo-FeNi, we adopted the 
primitive tetragonal structures represented by the solid 
lines in Fig. 1. The lattice parameter c along the per- 
pendicular (z) direction for each in-plane lattice constant 
a was optimized to minimize the total energy within a 
tolerance of 10 fieV. We checked that the MCA energy 
dependence on the in-plane lattice constants obtained us- 
ing the VASP-PAW code were well reproduced by all- 
electron calculations using the full-potential linearized 
augmentcd-plane-wave (FLAPW) method encoded in the 
WIEN2K packaged 



where S and L are the single-electron spin and angular- 
momentum operators, respectively. By expanding the 
\kncr) with an orthogonal basis of atomic orbitals labeled 
as /i(A), i.e., \kna) = J2iu c i^crNM (7 )j wc can obtain the 
second order contribution of i?so to the total energy as 
a sum over terms depending on spin-transition processes, 
atomic orbitals, and atomic sites, 
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where (a, <r']i,j) is an integral of joint local density 
of states (LDOS) given by, 
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where E-p is the Fermi energy. The matrix elements of 
(// t \L ■ S\/i t) arG given in Ref. HH as functions of the 
spherical angular coordinates 9 and 0, which are defined 
as polar and azimuthal angles between the spin quanti- 
zation axis and the crystal z axis. The MCA energy is 
derived as a difference in E^ 1 between in-plane (9 = n/2, 
<t> = 0) and perpendicular (9 = 0, </> = 0) magnetization 
directions, including the matrix elements (// j" \Lx |a* t) 
and (// t T)i respectively. We introduce a second- 

order contribution to the MCA energy depending on the 
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FIG. 2: (Color online) Bar graph of the second-order per- 
turbative contribution to the MCA energy depending on the 
atomic site and the spin-transition process of Llo-FeNi for an 
equilibrium lattice constant a = an = 3.556 A. 



atomic site and the spin-transition process as follows 
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We directly estimated the second-order perturbative con- 
tribution to the spin and atomic site dependent MCA 
energy , by calculating Cyj* using first-principles 

calculations without the spin-orbit interaction. c^ a can 
be obtained in the PAW formulation as a sum of the 
plane- wave part and the augmented part. Note that 
AE^p, in eq. (6) depends quantitatively on the choice 
of the spin-orbit coupling constant £$. We checked that 
the typical values £f g = 50 meV and £m = 100 meV— in 
eqs. (3) and (6) were consistent with the MCA energy 
and orbital magnetic moment of Llo-FeNi obtained from 
first-principles calculations that included the spin-orbit 
interaction. 



III. RESULTS 

A. Magneto-crystalline anisotropy and orbital 
magnetic moment 

We found that the magnetization of Llo-FeNi prefers 
the [001] direction, with an MCA energy of 0.078 meV 
per formula unit (f.u.) for the equilibrium lattice con- 
stants ao = 3.556 A and cq = 3.584 A, leading to K u = 
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FIG. 3: (Color online) Bar graph of the second-order per- 
turbative contribution to the MCA energy of nonvanish- 
ing angular momentum matrix elements between d states 

±\Wx W 2 G^' (a, I; i, i) and ±| <jj.'\L z \fj) \ 2 G^' (<T, U i, i) 
of Llo-FeNi for an equilibrium lattice constant a = ao = 3.556 
A. and"t=H" indicate spin-conservation and spin-flip 

terms, respectively. 



0.56 MJ/m 3 . This value is comparable with experimen- 
tal observations; i.e., 0.63 MJ/m 3 and 0.70 MJ/m 3 £i 
Note, however, that the long-range chemical order pa- 
rameters have been 0.5 or less for samples fabricated so 
far 4 On the other hand, for completely Llo-ordcred FeNi, 
the calculated order parameter value is 1. Experimental 
observations have indicated that the K u of Llo-FeNi is 
roughly proportional to the order parameter^ For Llo- 
ordered FePt, it was observed experimentally that K u in- 
creases as a quadratic function of the order parameter^ 
The same behavior was confirmed in previous theoretical 
studies based on the coherent-potential approximation 
(CPA)i 15 i 17 ' 20 The spin magnetic moments of the con- 
stituent atoms in Llo-FeNi are independent of the mag- 
netization direction, and were evaluated as 2.65 /xb and 
0.61 /xb for Fe and Ni atoms, respectively. On the other 
hand, the orbital magnetic moment of each atom shows 
a characteristic dependence on the magnetization direc- 
tion. The anisotropy of the orbital magnetic moment 
defined by A/i or b = fJ>^ — Mot™' takes a positive value, 
0.006 fiB, for Fe, whereas it is negative, -0.002 fiB, for Ni 
in Llo-FeNi. 

Recently, Kotsugi and his coworkers performed XMCD 
measurements of alternately layered FeNi thin films 
grown on CU3AU buffer layers^ They found that the 
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FIG. 4: (Color online) Local density of states (LDOS) of the d 
orbitals of Fe and Ni in Llo- FeNi for (a) a/ao = 1.00 and (b) 
a I ao = 0.904 without spin-orbit interaction as a function of 
energy relative to the Fermi energy. The equilibrium lattice 
constant ao is 3.556 A. 



angular distribution in the Fe orbital magnetic moment 
shows a clear and strong angular dependence, while the 
Ni orbital magnetic moment displays no clear angular 
distribution. This indicates that the anisotropy of the Fe 
orbital magnetic moment is a main contributing factor to 
the perpendicular MCA of Llo-FeNi, which is reasonably 
consistent with the present theoretical results. 



B. Origin of the perpendicular magneto-crystalline 
anisotropy 

According to Bruno's relation^ the MCA energy is 
proportional to the anisotropy of the orbital magnetic 
moment. To confirm a validity of Bruno's relation in 
Llo-FeNi, we evaluate the second-order perturbative con- 
tribution to the MCA energy AE^^ a , given in eq. (6). 
Figure 2 shows a bar graph of AE l a ^ a , for Llo-FeNi in the 
second-order perturbation. We found that the Fe atom 
makes a positive contribution to the perpendicular MCA, 
while the Ni atom makes a negative contribution. These 
results qualitatively agree with the behavior of the orbital 
magnetic moment. However, the spin-flip term AE^^ 
also makes a considerable contribution to the perpendic- 
ular MCA as does the spin-conservation term AEj^^. 
This means that the evaluation of MCA from the Fe or- 
bital magnetic moment using Bruno's relation^ will un- 
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FIG. 5: (Color online) (a) Total energy relative to that of 
a/ao = 1-00 (filled red diamonds) and the axial ratio c/a 
(filled blue squares), (b) the uniaxial MCA constant K u (filled 
black circles), and (c) the anisotropy of the orbital magnetic 

moments A/i or b = Morb 1 ' — Mort) ' °^ constituent Fe and Ni 
atoms calculated as a function of the in-plane lattice constant 
a for Llo-FeNi. The equilibrium lattice constant ao is 3.556 

A. 



derestimate the MCA energy of Llo-FeNi because of the 
absence of the spin-flip term in the orbital magnetic mo- 
ment. The large spin- flip term AE^ can be attributed 
to the weak exchange splitting of Fe between occupied 
majority-spin and unoccupied minority-spin states due 
to hybridization with Ni, resulting in a large majority- 
spin LDOS around Ep (see Fig. 4(a)). The small con- 
tributions of AE^ f and AEf^ t to the MCA energy 
can be attributed to the small LDOS of the unoccupied 
majority-spin states, which is typical properties of alloys 
including more-than-half transition metal elements. 

To obtain a further understanding of the origin of the 
perpendicular MCA, we show in Fig. 3 the second- 



order pcrturbative contribution to the total energy of 
nonvanishing angular momentum matrix elements be- 
tween d states, i.e., ±\(n'\Lx\n)\ 2 G^ (er, I; i, i) and 

±\(lJ,'\L z \n)\ 2 Gfy' (a,l;i,i). We found large matrix el- 
ements of the spin-conservation term between Fe d(xy) 
and Fc d(x 2 — y 2 ), which is the main contributing factor 
to the perpendicular MCA of Llo-FeNi. Furthermore, 
the spin- flip term from the majority-spin Fe d(3z 2 — r 2 ) 
to the minority-spin Fe d(yz) also makes a significant 
positive contribution to the perpendicular MCA. These 
results can be confirmed in the Fe LDOS of Llo-FeNi at 
a/ao = 1.00 in Fig. 4(a), where there are large LDOS of 
Fe d(x 2 — y 2 ), d(3z 2 — r 2 ), and d(yz) in the vicinity of 



C. Effect of tetragonal distortion 

We examined the dependence of A^ u on the in-plane lat- 
tice constant a of Llo-FeNi, in which the perpendicular 
lattice parameter c is relaxed for each in-plane lattice con- 
stant a. The change in the axial ratio c/a and the total 
energy are shown in Fig. 5(a). As shown in Fig. 5(b), K u 
increases with decreasing in-plane lattice constant (with 
increasing perpendicular lattice parameter c), reaching 
a maximum value of 1.6 MJ/m 3 at a = 3.182 A. (c/a 
= 1.242). The present result suggests that the compres- 
sive in-plane stress utilizing the lattice mismatch between 
Llo-FeNi and the buffer layers is effective in achieving a 
higher MCA energy. Indeed, previous experiments have 
indicated that K u monotonically increases with increas- 
ing axial ratio c/a of Llo-FeNi^ According to previous 
first-principles calculations,— the MCA energy of Llo- 
FeNi reaches a maximum value of 0.90 MJ/m 3 at c/a 
= 0.95. The discrepancy between this result and that 
obtained in the present study may partly be attributed 
to the computational method adopted in the previous 
study, i.e., the linear muffin-tin orbital method within 
the atomic sphere approximation, which assumes a spher- 
ically symmetric potential in each atomic sphere. 

To elucidate these results, we show in Fig. 5(c) the 
anisotropy of the orbital magnetic moments A/i^ h and 
A/!^ between cases in which the magnetization is along 
the [001] or [100] direction as a function of in-plane lat- 
tice constant a. A/j,^ h decreases with decreasing in-plane 
lattice constant from a/ao = 1-00, then slightly increases 
with increasing a/ao. On the other hand, A^^! b de- 
creases with decreasing in-plane lattice constant from 
a/ao = 1-00. These behaviors are inconsistent with the 
dependence of the MCA energy, indicating that the spin- 
conservation term is not a primary contributing factor 
to the large perpendicular MCA of Llo-FeNi under com- 
pressive in-plane distortion. To understand this point, 
we show in Fig. 6 the spin-resolved MCA energy AE^^ 
and AE'^^ as a function of the in-plane lattice con- 
stant. First, the spin-conservation terms AEj^ and 
A-E^i, reproduce the in-plane lattice dependence of the 
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FIG. 6: (Color online) The second-order perturbative contri- 
bution to the MCA energy AEn _j, (i=Fe or Ni) as a function 
of in-plane lattice constant a for Llo-FeNi. The equilibrium 
lattice constant ao is 3.556 A. 

difference in orbital magnetic moments A^ h and A^ h 
shown in Fig. 5 (c), because the orbital magnetic moment 
originates only from the spin-conservation term. Further- 
more, the spin-flip terms AE^^ and A£^A^ increase 
with decreasing in-plane lattice constant, which is con- 
sistent with the dependence of the MCA energy shown in 
Fig. 5(b). This means that the large perpendicular MCA 
energy of Llo-FeNi under compressive in-plane distortion 
can be attributed to an increase in the spin-flip terms of 
Fe and Ni (AE^ i and Afi^). Figures 7(a)-(h) show 
the MCA energy of nonvanishing angular momentum ma- 
trix elements between d states ±|(/x'|Lj> ( :|A t )| 2< -'^ A [ 1 (er, 4, 
and ±\([i'\L z \/j,)\ 2 G l / ll l [ l (a,l;i,i) as a function of 
the in-plane lattice constant. It was found that the posi- 
tive contribution to the MCA energy decreases, while the 
negative contribution increases, with increasing in-planc 
lattice constant. Although these contributions almost 
cancel each other out, there is a net increase in MCA 
energy with decreasing in-plane lattice constant. Among 
them, the matrix elements with (a; 2 — y 2 \L Z (x)\xy) and 
(xy\Lz(x)\ x2 ~ y 2 ) strongly depend on the in-plane lat- 
tice constant, in both the spin-conservation and spin-flip 
terms. We conclude that the significant reduction of the 
negative contribution related to the matrix elements of 
(x 2 - y 2 \L z(x) \xy) and (xy\L z{X ) W ~ y 2 ) lead to the 
increase in the MCA energy of Llo-FeNi with decreasing 
in-plane lattice constant. We confirmed that the LDOS of 
d(x 2 — y 2 ) and d(xy) are remarkably modulated by com- 
pressive in-planc distortion compared with the LDOS of 
other orbital components. As can be seen in Figs. 4(a) 
and (b), the LDOS of the majority-spin d(x 2 — y 2 ) below 
E-p are reduced by the tetragonal distortion from a/ao 
= 1.00 to a/ao = 0.905. This can be attributed to an 
enhancement of the delocalization of the d(x 2 — y 2 ) and 
d(xy) orbitals resulting from strong bonding between the 
in-planc d orbitals under the compressive in-plane distor- 
tion. 
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FIG. 7: (Color online) The 3d orbital contribution to the 
MCA energies of Llo-FeNi as a function of the in-plane lat- 
tice constant a for Fe (a)~(d) and Ni (e)~(g). Here, we 
use abbreviate notations for the MCA contribution of the 
nonvanishing angular momentum matrix elements between 
d states, i.e., ±(/x' 4- \Lz\n&) (" + " corresponds to a =4- 
and " - " to a =f) for ±| {n'\Lz |m> \ 2 G^'/ (a, I; i, i), and 
±(fi' I \L x \v&) (" + " corresponds to a =t and " - " to 
a —\.) for ±|(^i'|Lx|/i)| 2 Gu / / t (a, The equilibrium lat- 

tice constant ao is 3.556A. 



IV. DISCUSSION 

A. Effect of orbital-polarization 

It is well known that the orbital polarization energy 
related to the second Hund's law, which is not consid- 
ered in the present calculation, enhances the MCA ener- 
gies of most metals and alloys. According to Ravindran 
et al.^ the MCA energy of Llo-FeNi is enhanced from 
0.077 mcV/f.u. to 0.172 meV/f.u. due to the orbital- 



polarization energy. To clarify the effect of the orbital 
polarization (OP) on the MCA energy in the tetragonal 
distortion of Llo-FeNi, we performed calculations of the 
MCA energy including the orbital polarization energy us- 
ing the FLAPW method with the WIEN2k code£I using 
the GGA-PBE function.^ Figure 8 plots the MCA en- 
ergy and the anisotropy of the orbital magnetic moment 
with and without OP as a function of the in-plane lat- 
tice constant. At the equilibrium lattice constant (a/an 
= 1.00), we found that the MCA energy with OP is 
roughly two times that without OP, which is consistent 
with previous calculation resultsJ^ The origin of the en- 
hancement was well discussed in Ref. [13, i.e., the or- 
bital polarization energy effectively enhances the spin- 
orbit coupling parameter = + BL, where B is a 
Racah parameter— expressed in terms of Slater integrals 
of the single-particle wave functions for the d orbitals. 
Furthermore, we found that the MCA energy with OP in- 
creases with decreasing in-plane lattice constant, which is 



(a) 




3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 
In-plane lattice constant a [A] 




3.2 3.3 3.4 3.5 3.6 3.7 3.! 
In-plane lattice constant a [A] 



3.9 



FIG. 8: (Color online) (a) The uniaxial magneto-crystalline 
anisotropy (MCA) energy K u and (b) the anisotropy of the 
orbital magnetic moments of constituent Fe and Ni atoms 

A/U rb = /vt> — Morb°' with and without the orbital polar- 
ization energy as a function of the in-plane lattice parameter 
a /ao calculated by the FLAPW method^. The equilibrium 
lattice constant ao is 3.556 A. 
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similar to the results without OP. However, the enhance- 
ment of the MCA energy under the compressive in-plane 
distortion is suppressed compared with the MCA energy 
without OP. The MCA energies with OP around a/ao = 
0.90 are almost identical to the MCA energies without 
OP, despite the increase in the anisotropy of the orbital 
magnetic moment caused by the orbital polarization en- 
ergy, as shown in Fig. 8(b). This can be better under- 
stood by considering the different origin of the perpen- 
dicular MCA of Ll -FeNi under tetragonal distortion. 
As was discussed in the previous section and shown in 
Fig. 6, the perpendicular MCA around the equilibrium 
lattice constant (a/a = 1.00) is mainly caused by the 
spin-conservation term of Fe (ASj^), while the spin-flip 
terms of Fe (A£&. , ) and Ni (Ai^j^) provide the largest 
contribution under highly compressive in-plane distortion 
(a/a <0.95). The orbital polarization energy increases 
the orbital magnetic moment and its anisotropy, which 
is equivalent to an enhancement of the spin-conservation 
term, but it does not enhance the spin-flip term because 
of the absence of a spin operator. Therefore, the effect of 
the orbital polarization energy is remarkable in the MCA 
energy arising from the spin-conservation term compared 
to that from the spin-flip term. 

B. Transport properties through the MgO barrier 

Finally, we discuss the potential of Llo-FeNi as a fer- 
romagnetic electrode for magnetic tunnel junctions with 
an MgO barrier. Figure 9 shows the electronic band- 
structures along the [001] direction of bulk Llo-FeNi. The 
totally symmetric band in FcNi crosses the Fermi level in 
both the majority-spin and minority-spin states. The 
minority-spin Ai band is predominantly composed of Fe 
and Ni d(3z 2 — r 2 ) orbitals at E-p, while the majority- 
spin band is mainly constructed from Fe and Ni p(z) 
orbitals. Similar results have been obtained for the 



(a) Majority-spin state (b) Minority-spin state 




0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 



k z at (/t„/t y )=(0,0) [2it/a] K at (k x ,k y )=(0,0) [2n/a] 

FIG. 9: (Color online) Electronic band-dispersion curves rel- 
ative to the Fermi energy along the [001] direction of bulk 
Llo-FeNi for (a) majority- and (b) minority-spin states. 



band dispersion of Llo-FePti 3 - 7 - Since the Ai band elec- 
trons predominantly transmit the MgO barrier ; 38 i 39 huge 
TMR ratios cannot be expected in Llrj-FeNi/MgO/Llo- 
FcNi(OOl) MTJs from the viewpoint of symmetry com- 
patibility between the electronic band structures in bulk 
FeNi and MgO. In fact, we obtained a very small tunnel- 
ing magnctorcsistance ratio of 51% for the Fe-terminatcd 
Ll -FcNi/MgO/Ll -FeNi(001) MTJ (11% for the Ni- 
tcrminatcd MTJ) using Landauer-type ballistic conduc- 
tance calculations. For details regarding these conduc- 
tance calculations, see ref. H3 and the references therein. 
Thus, to achieve higher TMR ratios, it is necessary to in- 
sert appropriate ferromagnetic layers, e.g. Fe or CoFeB, 
between the Llo-FeNi electrodes and the MgO barrier 4£ 



V. SUMMARY 

We evaluated the MCA energy of Llo-FeNi using first- 
principles density-functional calculations. The perpen- 
dicular MCA found in Llo-FeNi can be attributed pre- 
dominantly to the constituent Fe atoms. The MCA en- 
ergy of Llo-FeNi increases with decreasing in-plane lat- 
tice constant. The perpendicular MCA under tetragonal 
distortion was elucidated as follows. The perpendicular 
MCA at the equilibrium lattice constant is mainly caused 
by the spin-conservation term in the second-order pertur- 
bation of the spin-orbit interaction. On the other hand, 
the perpendicular MCA under highly compressive in- 
plane distortion can be attributed to the spin-flip terms 
of Fe and Ni. The orbital polarization energy increases 
the MCA energy with a small tetragonal distortion due 
to enhancement of the spin-conservation term, while the 
effect of the orbital polarization energy is suppressed un- 
der highly compressive in-plane strains. We concluded 
that modulation of the in-plane lattice parameter of an 
Llo-FeNi thin film by choosing appropriate buffer layers 
is effective in obtaining large MCA energies correspond- 
ing to 1.0 MJ/m 3 , and is worth further investigation. 
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